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Abstract 

We consider a p-wave superconductor wire coupled to two metallic rings, confined to a one- 
dimensional wire. At the two interface between the the wire and the metallic rings the pairing 
order parameter vanishes, as result two zero modes Majorana fermion appear. The two metallic 
rings are pierced by external magnetic fluxes. The special features of the Majorana Fermions can 
be deduced from the correlation between the currents in the two rings. 
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I Introduction 



Topological Superconductors are characterized by the invariance under charge conjuga- 
tion symmetry. As a result of this invariance zero modes Majorana fermion appear at the 
interfaces between the superconductor a metal. Sr 2 Ru is the material where Majorana 
fermion might be observed since the pairing order parameter is characterized by p x + ip y 
symmetry. The electronic excitations for the ground state pairing p x + ip y are given by half 
vortices which are zero mode Majorana fermions 

The realization of the p-wave superconductors physics and the formation of the zero mode 
Majorana Fermions at the edges of the wire can be achieved using a p-wave wire of length 
L coupled to two metallic rings which are pierced by external magnetic fields. The current 
in the rings is coupled to the p-wave wire trough the Majorana modes which are bound 
at the interface. For different fluxes in the rings we find that the excitations in the wires 
become imaginary resulting in an unstable vanishing current. The only stable current are 
obtained for the case that the imaginary part vanishes. The stable current are obtained for 
special relations between the magnetic fluxes and the wire excitations energy. This feature is 
attributed to the existence of the Majorana fermions. We find that the correlation between 
the currents in the two rings is affected by the presence of the Majorana Fermions. Therefore 
the model introduced here can be used for the identification of the Majorana fermions. 

The content of this paper is as follows: In chapter II we present the model of the the 
p — wave wire coupled to two metallic rings. Using the left and right mover we obtain 
the continuum representation of the superconductor wire. At the two edges of the wire we 
obtain the two zero modes of the wire. The effect of the coupling between the wire and the 
metallic rings is dominated by the zero modes of the wire. As a result we derive the effective 
Hamiltonian between the zero modes and the metallic rings. This effective Hamiltonian is 
controlled by the low energy excitation in the wire (the coupling energy between the two 
edges) e ~ |Ao|e~ L ' A °' ,L|Ao| >> 1 where Ao is the pairing field and L is the length of the 
wire. In section /// we consider the case that the wire energy e — > 0. In section IV we 
consider the case where the wire energy is finite. Section V is devoted to conclusions. 

II- The model for the p-wave wire weakly coupled to two rings pierced by 
external fluxes 
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The p-wave wire of length L is given by the Hamiltonian Hp-w '■ 
Hp_ w = -tJ2(C + (x)C(x')+h.c.)-^ F J2C + (x)CXx)-AJ2h^ 

(1) 

The pairing gap is given by A and the polarized fermion operator is given by C(x) = C a= ^{x). 
The matrix elements 7 XjX / obey obey the p-wave symmetry: j x ,x'=x-a = —lx,x'=x+a, \lx,x'\ — 
1, therefore the time reversal and parity symmetry are both broken and obeys the pairing 
boundary conditions A(x = 0) — A(x — L) = and A(x) = A for < x < L. 

We introduce the right and left fermions in the continuum representation for the fermions 
in the wire : C(x — na)y/a — > C(x) = e lkpX C R (x) + e~ lkFX C L (x) and find that equation 1 is 
replaced by the Hamiltonian : 

H P _ W = J dx[v F ^(x)a z {-id x )^(x) + A(x)^(x)a x ^{x)} (2) 

This Hamiltonian is invariant under the charge conjugation symmetry. The spinor ^(x) = 
[C r (x),Cl(x)] = [e~C R (x), c~Cl{x)] satisfies the reality constraint condition K^(x) = 
ty(x), where K is the charge conjugation operator. 

The pairing field A{x) = AA(x) sm(kpa) can be written as A(x) = Ml(x) + Mp(x — L) 
where M L (x) = ^p-sgn(x) and Mp(x-L) = ^p-sgn(x-L) obeys the domain wall property: 
M L (-x) = -M L (x) (at x=0) , M R (-(x - L)) = -M R (x - L) (at x=L) . The zero modes 
eigenfunctions are given by r}\(x) = [r]i(x), r] 2 (x)] T and eigenstates of the operator a y rj\(x) = 
\r]x(x) with A = ±1. The zero mode spinor which is localized around x = is identified 
with A = — 1 and the second one which localized around x = L is identified with A = 1 

VLeft(x) = 77 A= _i(a;) = Jo [ X> _ t y 



i \ i \ — f?Mx')dx'e l 4 



(3) 



The spinor operator ty(x) with the two zero mode Majorana operators ai (at the left edge) 
and a r (at the right edge), (a r ) 2 = (a/) 2 = \ takes the form : 

->■ V(x) + a r r] X=1 {x) + a^-iix) (4) 

As a result the low energy of the p-wave wire is given by: 

Hp_ w = J dx[v F ¥(x)a z (-id x )^(x) + A(x)¥(x)a x ^(x)] « % -eaia r (5) 
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where e « |A |e~ L l Ao l ,L|A | >> 1. 

At this stage we include the two rings Hamiltonian pierced by the fluxes <fi , i — 1, 2 and 
length / r j„ 9 « L. The left ring is restricted to the region —l r i ng < x < and the right 
ring is restricted to L < x < L + l r i ng . Since only the wire fields at x = and x = L are 
involved we fold the space of the right ring % = 2 such that both rings are restricted to the 
region —l r ing < x < 0. As a result the external fluxes obey: ipi — > (pi and ip 2 — > — y?2- In 
addition we replace for each ring the fermion operator ^ (a?) ,i=l ,2 by the rig/it -R«(x) and 
left fermions Li(x): 

il>i(x) = Ri(x)e ikFX + U[x)e- lkFX (6) 
We replace for each ring the right and left movers by four Mayorana operators: 

ri = (Ri(p)-I$(p))(-i) 
li = Li(p) + L\(p) 

(7) 

The the matrix element between the wire and rings is given by —g. As a result the low 
energy hopping Hamiltonian is given by: 

H T = ^[a l (r 1 + l 1 ) + a r (r 2 -l 2 )] (8) 

We replace the Majorana zero modes «/ and a r by the fermion pair, q = ai+ia r , q^ = ai—ia r 
which obey: [q, q^] + = 1, q^\0 >= |1 > and q\l >= |0 > where |0 > is ground state of wire 
and rings: Ri, p {x)\0 >= L iyP {x)\0 >= R[ h {x)\0 >= L[ h (x)\0 >= q\0 >= 0. Where Ri, p {x), 
Li, p (x) are the particle operators and R ih (x), L\ h (x) are the holes operators R\ h (x), L\ h (x). 
The right and left mover are given as a linear combination of particles and holes operators. 
Ri tP (x) L i>p (x) represent the annihilation of particles and R\ h (x), L\ h (x) are the creation 
operators for holes. 

R t (x) = R hP (x) + R[ h (x); L t (x) = L hP (x) + L\ h {x) (9) 

Using the Fermionic representation we replace Hp_w given in equation 5 and Ht given in 
equation (8) by : 

H P _ W + H T = eq^q - ^=[(q + q^( n + h) - i(q - J)(r 2 - l 2 )\ (10) 
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The value of the wire energy e in equations 5, 10 is based on the projection of the spinor (in 
equation 4) on the zero modes rjx=i(x) and n\=-x(x). Since the leads couple to the modes 
in the wire, we expect that the non-zero modes will give rise to a finite width of the energy 
e. Therefore for finite energies we will replace e by e = e — iT, where the width T oc g 4 . 

We perform an exact integration over the fermion operators q,q^ and find the time de- 
pendent effective interaction H e ff(t): 

Heffit) = ^f J dtfi[t,1f\e*lW[r 3 (t) - k(t) + i(n(t) + li(t))][r a (f ) - W) - i(n(t') + h(t'))] 

(11) 

where //[£,£'] is the step function which is one for t > t' . 

III-The effective interaction H e ff(t) in the limit e — > 

When e — > equation 11 is replaced by : 



Heff(t) 



-w 



dt'fi[t,t'}[r 2 (t) - l 2 {t) + i{ Tl {t)) + h(t)][r 2 (f) - h{t') - i(n(0 + )] 

(12) 



Using the scaling analysis given in 12[ we observe that the effective interaction flows to 



the strong coupling limit and find g 2 (b) = g 2 b 



2L2-Q 



a 



1 where b > 1. Since the coupling 



constant <?(&) flows to infinity, the only way a solution will exists if the effective interaction 
annihilates the ground state: H e ff(t)\0 >= 0. Therefore the physical solution is given by 
the constraint |4j equation: 



fa® - h{t) - i(n{t) + h{t)]\0 >= 



(13) 



Since Ri iP (x)\0 >= L iiP (x)\0 >= R\ h (x)\0 >= L\ h (x)\0 >= the constraint condition 
implies for particles (p stands for particles) the equation : — M. P ) ~~ (^Ip — Ri p )]\® >— 

0; and for holes (h stands for holes) [i(R2,h ~ Li t h) — (L 2t h — Ri,h)]\0 >= 0. We find the 
constraint equation: 



^(x = 0) = R x {x = 0)+L 1 (x = 0) = e~^[R 2 {x = 0)+L 2 (x = 0)] = e~^ijj 2 {x = 0) = ^ 2 {x = 0) 

(14) 
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The explicit identity contains the phase factor e ! 2 which is obtained from Bosonization 
(see below). Following equation 14 we find the constraint condition for the ground state 
|0 >: 

k= [ipi(x) - ip 2 (x)]\ x=0 , k\0>=0 (15) 
Following 4] we find two additional constraints equation : 

£= \(-id x - ^^fM*) ~ Hd* + ^02) 2 Mx)]\x=o, £\0>=0 (16) 

''ring ''ring 

2tt 2tv 
J=[(-id x -- 1 )^ 1 ( x ) + (-id x + - 2 )^ 2 {x)\ U=o, J|0>=0 (17) 

''ring ''ring 

Any eigenstate of N particles must satisfy the set of equations 15 — 17 with the periodic 
boundary condition ifji(x) = vf)i(x + Uing)- For the case N = 1 (one particle) we have : 

\N = 1>= f° dx[h{x)4{x) + f 2 (x)4(x)} (18) 

J Iring 

where fi(x) are the amplitudes . Using equations 15 — 17 we find finite solutions for the 
amplitudes fi(x) only when the fluxes are equal. 

IV-The finite limit e ^ 

In order to study the finite limit e ^ we will use the zero mode Bosonization method 
[si Q]- The right Ri(x) and left Li{x) fermions for each ring i — 1, 2 is given by: 

r.( x ) = ./A^e^e^^'^^^W 



2vr 

(19) 



Li{x) = ^A^e^^W^-' % x e iV*irtLA*) 



Where ZiZ 2 = —Z 2 Z\ are Majorana variable which ensure the anti-commutation between 
the two rings in the bosonic representation. In the bosonic representation we have the 
zero modes a R j, bosons and their conjugates ,N R j,N L j. The zero modes obey the 
commutation rules :[— a^, Nlj] = iS^j and [a Rt i, Nrj] = iSij. 

As a result we obtain the zero mode representation in terms of the fermion numbers 
Ni,i,Nji t i and fluxes (pi in each ring: 

#o = ^[(A^i - N R>1 + 2^) 2 + (N L}1 + N RA ) 2 ] + ^{(N L>2 - N R>2 + 2<p 2 f + (N Lfl + N R>[ 

^ring ^ring 



(20) 



Using the equations of motion ifo-^f 1 = [oiR t i, H ] and ifo—^- = [&L,i, Ho] ,i — 1, 2 we obtain 
the zero mode representation in the interaction picture. We will use the zero mode 
representation a z Ri (t), a T Li (t) in the interaction picture in order to evaluate equation 15. 
We find that H e ff(t) is given in terms of the zero mode functions F(t + |) and G{t — |) : 

drF(t+-)e-^G(t--) (21) 

where g 2 = g 2 ^- is the coupling constant and u = f , z/j = j^Ai are the equivalent wire 

Z7T /I ''ring 

and rings frequencies and f = ^ is the width. The functions F(t + |) and G(t— |) are given 
by: 

F(t + ^) = Z 2 ((sin4 j2 (t + I) - cosai j2 (t + ^)) + iZ 1 ((sm a^t + I) - coso^f + ^)) = 
Z 2 [sin a^ 2 (t) cos(z/ 2 r) + COSQ 4, 2 (^) sin(z/ 2 r) — cosa^ 2 (t) cos(z/ 2 r) + sino;^ 2 (t) sin(^ 2 r)] 
+i.Zi[sin aR i(t) cos(^ir) — cosa^,-^) sin^r) + cosoj£ 1 (t) cos(z/it) — sina^^t) sin^r)] 

(22) 

G (* - = Z 2((sina^ 2 (t - -) - cosa£ 2 (t - -)) - ^((sine^t - -) + cosoi^t - -)) = 
Z 2 [sin (*) cos(z/ 2 r) — coscr^ 2 (t) sin(z/ 2 r) — cosa£ 2 (t) cos(z/ 2 r) — sina£ 2 (t) sin(z/ 2 r)] 
— zZi[sina^ 1 (t) cos(z/it) + cosa£. 1 (t) sin(z/ir) + cosa£ 1 (t) cos(z/it) — sina^^i) sin(^ir)] 

(23) 

We perform the integration with respect r and find: 

H eff « fl^ + (24) 
where * s the rea ^ P ar ^ °f the effective action: 

zjreal 
H eff = 

' : [cos(2a fi)2 ) - cos(2a Lj2 ) - cos(2a Lj i) + cos(2a fl> i) + 2(sin(a fl>2 + «l,2) - sin(a fl> i + c*l,i))] 



^o 2 + r 2 



(z/ - /a) 2 + r 2 (i/ + ^i) 2 + r 2 



(i/ - ^ 2 ) 2 + r 2 (z/ + z/ 2 ) 2 + r 2 



(25) 
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The imaginary part H^fj of the action causes the current to vanish. Finite solutions will be 
obtained for the ground states |0 > which obey H^\0 >= 0. Therefore the solutions for 
finite currents are equivalent to a constraint condition for the ground state |0 >. 



r 



r 



H eTf = [- — 

[(sin(a i?j2 ) - cos(a Lj2 ))(cos(a Li i) + sm(a Rjl )) + (sin(a flj2 ) + sin(a £|2 )) 



2 / 



r 2 



(z/ + ^=^ 2 



r 2 



(cos(a Rj i) + sin(o!L,i))] 
+Z 2 Z^ 2 [- 



r 



> -^) 2 + r 2 (z/ + ^) 2 + r 2J 

[(sin(o!fl ) 2) - cos(a Lj2 ))(cos(a L ,i) + sm(a Ril )) - (cos(a jRj2 ) + sin(a Lj2 )) 
(cos(a R) i) - sin(a Ljl ))] 



We find two cases for which a finite solution exist. 

The first case corresponds to v (\ u ' 2 ^ Ul |) with the solution: 

Hi%\0>=0 

«R,2 = OR,l + U 

The second case corresponds to u ~ ( i ^jr 1 ) with the solution: 

^/7I°>= 

«i?,2 = -Q!fl,l + U 
«L,2 = -Q!L,1 - U 

We introduce the definitions for the zero mode fields : 

an = a 2 = a; fa — fa — 2u, fa = ft 



(26) 



(27) 



(28) 



(29) 



The solutions are independent on the arbitrary field u which plays the role of a gauge 
condition and has to be integrated out. We integrated over the field u we find the effective 
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Hamiltonian for the conditions z/ ss ( U2 2 Ul ) ana ^ u o ~ ( l/2 1 2 Ul )- We introduce the 
magnetic flux the variables (p and A: 

V= g ' 2 ^ ^ 

We can write both cases in a closed form: 

H irvp d A , 9 . <i , 9l 7rvp r , d . , , d , 9 , 
¥ ^ 2l- [( -^ + ^ + ( -^ } 1 + 2l- [( "^ + ^ + ( -^ } 1 + 

r sin( a )sin(^)-2^ h + 7 , + . ^ , J ^ 



L z/ 2 + r 2 > -(<^ + A) 2 + r 2 (i/ + (^ + A) 2 + r 2j 

o~2r ^Q-(^ + A) Z/q + (g + A) Z/q - (g - A) Z/q + (g - A) 

>o-(y? + A) 2 + f 2 (i/ + (<p + A) 2 + f 2 (u - (tp - A) 2 + f 2 (i/ + — A) 2 + f 2 



(31) 



The first line of equation 31 represents the the Hamiltonian for the two metallic rings pierced 
by the external fluxes expressed in therms of the zero mode of the metallic rings. The second 
part of equation 31 represents the coupling between the wire and the two rings. We observe 
that this part is restricted by the constraint condition ifl ~' P2 = i/ or y ' 1 "^ 2 = vq. This 
constraints represents the effect of the Majorana fermions on the p-wave wire. 

In order two investigate the Hamiltonian in equation 31 we will use the algebra of the 
zero modes {J [3, Q: 

J = N R -N L = ~2i^-; Q = N R + N L = -2z^ (32) 
da dp 

with the eigenvalues and commutation rules: 

J\J, Q >= J\J,Q >; J = 0, ±1, ±2, .. 
Q\J, Q>=Q\J,Q>;Q = 0, ±1, ±2, .. 

From the commutation relations [a, J] = 2i, [/3, Q] = 2i we establish the relations : 

e ia \J,Q >= \ J+ 1,Q >;e~ ia \J,Q >= \J-1,Q> 
e ip \J,Q>= \ J,Q + l >;e~ if} \J,Q >= \ J,Q-1> 



(33) 



(34) 



9 



The eigenf unctions are given by: 



< a\ J, Q = >= ^=e iaJ ] < P\ J = 0, Q >= -Le* Qj 

V47T V47T 



(35) 

Using the algebra of the zero modes we compute to lowest order (in perturbation theory) 
the energy for the ground state of the two rings coupled to the wire. As a function of 
the coupling constant A = < 1 and maximum frequency v max which is given by the 
electronic bandwidth frequency. We find for the ground state energy E(ipi, (p 2 ): 

E{0i, 02) = + A) 2 + {<p — A) 2 ) 

''ring 

w *W w Vp - (<p + A) Up + (g + A) 

^ (z/ - ((^ + A) 2 + f 2 + (y? + A) 2 + f 2 



v -((p-A) up + (<p- A) 

+ 7 — rv^ — * \ o , aJ JCW + <Wa) 



(i/ - - A) 2 + r 2 (v + (<p- A) 2 + r 2 ' 



(36) 



Using eq.36 we compute the currents Ii = dE ^^ for the two rings i — 1,2 using the 
conditions: z/ ps £i_^2 = ^ 

The current in ring one I\ and ring two J 2 are represented in terms of (p — LplJ ^ LP2 ; z/ = A 
and the current amplitude I = ^p^- 

'•ring 



h x r ^ 2 -r 2 -(2z, + v? ) 2 + f 2 

— = </? + f — A - 1 - 

h V + r 2 ) 2 (2// + V 9) 2 + r 2 ) 2j 

/ 2 _ u p 2 -P -(2^ + ^) 2 + f 2 1 

— = ip — UQ — A 7 1 7 

A) V + r 2 ) 2 (2// + V 9) 2 + r 2 ) 2j 

^=[^-A( + - (2 ^ + ^ + f2 )] 

2/ V + r 2 ) 2 (2z/ + 0? + r 2 ) 2 



(37) 



In figure 1 we have plotted the current Jl 2 "^ 2 as a function of the flux Cp for the case ~ 
ipl ~ ip2 = A — 0.01. We observe that for (p > 0.1 the current is proportional to (p. From 
other-hand when the (p < 0.1 the current in each ring is affected by the flux in the other ring. 
This is seen from the negative contribution of the current shown in figure l.The negative 
current contribution might be related to the Andreev reflection which occurs at the interfaces 
between the superconductor and the metal. 
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In figure 2 we have plotted the current as a function of the flux (p for the case e = 
considered in chapter III. Due to the constraint condition is = yi ~ y2 we have the relation 
uq = Lpl ^ Lf " 2 = A = 0.. We find a stable current in agreement with chapter III where the 
current scales linearly with the flux (p\ — (p^ — <fi- 

V-Conclusion 



We have inv estig ate the dependence of the current on the fluxes for the entire regime of 
parameters 
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15] In the limit of large L , e — > the current vanishes in both rings when 



the two fluxes are different. We observe that for a finite energy e and different fluxes the 
current dependence is more complex. When the two fluxes are almost equal the current is a 
function of the averaged flux. For the case that the flux difference is comparable to the flux 
average, the current changes sign. We can interpret this effect as an Andreev reflection and 
represents a finger print of the Majorana fermions. 
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